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Abstract 

The first hitting time process of an inverse Gaussian process is considered. It is shown that 
this process is not Levy and has monotonically increasing continuous sample paths. The 
density functions of one-dimensional distributions of the process are obtained. Its distribution 
functions are not infinitely divisible and their tail probability decay exponentially. It is also 
shown that the hitting time process of a stable process with index < /3 < 1 is not infinitely 
divisible. The density function is shown to solve a fractional partial differential equation. 
Subordination of the hitting time process to Brownian motion is considered and the underlying 
PDE of the subordinated process is derived. 

Key words: Hitting times; infinite divisibility; inverse Gaussian process; tail probability, sub- 
ordinated process. 



1 Introduction 

The first hitting time process (first passage time or first exit time) is a stopping time process 
which arises naturally in diverse fields such as finance, insurance, process control and survival 
analysis (Ting Lee and Whitmore (2006)). When the total claim size S{t) of an insurance 
company follows inverse Gaussian process (see e.g. Dufresne and Gerber (1993); Gerber (1992) 
and references therein), then the hitting time process corresponds to the situation when the 
company's aggregate claim exceeds certain level. Indeed, this is closely related to the proba- 
bility of ruin also. As another example, the first time an asset price process reaches a certain 
level is of interest to an investment firm. Hitting-time models are also be applied to expected 
lifetimes of mechanical devices, where the device breaks down when the process reaches an 
adverse threshold state for the first time. 



For the standard Brownian motion the hitting time process has been well studied in the 
literature and the hitting time process called Levy subordinator (see e.g. Applebaum, 2009). 
Hitting time for a general Gaussian process is obtained by Decreusefond and Nualart (2008). 
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Many authors have discussed the hitting time process for strictly increasing Levy processes (see 
e.g. Meerschaert and Schcfilcr (2008), Veillette and Taqqu (2010a)) and references therein). 
However, an cxphcit expression for the density of hitting time process is not possible in many 
cases (see e.g. Veillette and Taqqu (2010b)). The focus of many of these authors are on the 
Laplace transform of hitting time density or on mean hitting time process. In this paper, our 
aim is to discuss the hitting times for an inverse Gaussian (IG) process, which has not been 
addressed in the literature and investigate its various other properties. 

The density function of an inverse Gaussian distribution IG(a, 6), with parameters a and b 
(a, 6 > 0), is given by 

fix; a, b) = (27r)-'/'ax-='/'e"''-^(«'^"'+'''^), x > 0. (LI) 

It is well known that inverse Gaussian distributions, and generalized inverse Gaussian (GIG) 
distributions, are generalized gamma convolutions (see Halgreen, 1979) which arc infinitely 
divisible. The IG process {G{t),t > 0} with parameters 5 > and 7 > is defined by (see 
Applebaum, 2009, p. 54) 

G{t) = inf{s > 0; B{s) + 7s > 5t}, (1.2) 

where {B{t),t > 0} is the standard Brownian motion. That is, G{t) denotes the first time 
when a Brownian motion with drift 7 hits the barrier 6t. The IG subordinator G{t) is a 
non- decreasing Levy process such that the increment G{t + s) —G{s) has an inverse Gaussian 
IG(5t,7) distribution. The Levy measure tt corresponding to the inverse Gaussian process is 
given by (see e.g. Jeanblanc et.al. (2009), p. 596) 

T^idx) = --j==e'~I{:c>o}dx. (1.3) 

Note that almost all sample paths of {G{t)} are strictly increasing with jumps, since the sample 
paths of {B{t) +7^} are continuous and having intervals where paths are decreasing. In other 
words, < t2 =^ ^(^1) < ^'(^2) almost surely. This can be seen as follows. Consider 
= {s > : -B(s) +75 > 5ti} and At^^ {s>0 : B{s)+js > 5t2}. This implies for ti < t2, 
^ and hence G{ti) < G{t2). Next suppose that G{t) is not strictly increasing, this 
means that for some ti < t2 we have G{ti) — G{t2). But due to continuity of sample paths 
of BM with drift, G{ti) = ^(^2) implies B{G{ti)) + -fG{ti) = = B{G{t2)) + lG{t2) = ^2 
and which is a contradiction. The strictly increasing property of sample paths of {G{t)} also 
follows by using Theorem 21.3 of Sato (1999), since 

/•oo 

/ 'K{dx) = 00. (1.4) 
Jo 
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Let H{t) denote the right continuous inverse of G{t), defined by 

H(t) = mi{u > : G{u) >t},t> 0. (1.5) 

Since the sample paths of G{t) are strictly increasing with jumps, the sample paths of H{t) 
are continuous and are constant over the intervals where G{t) have jumps. In this paper, we 
find out the Laplace transform of the density of the process H{t) with respect to the time 
variable and then invert it to get the corresponding density function h{x,t). It is shown that 
the process H{t) is not a Levy process. The distribution of H{t) is not infinitely divisible and 
its tail probability P{H{t) > x) decays exponentially. The density function of h{x,t) solves a 
pseudo fractional differential equation. The PDE satisfied by the subordination of H[t) with 
standard Brownian motion B{t) is also obtained. 

2 Hitting Time of Inverse Gaussian Process 

Let Ct'w{x,t) = e~'^^w{x,t)dt denotes the Laplace transform (LT) of the function w with 
respect the the time variable t. For a Levy process D[u) with corresponding Levy measure ttd 
and density function / from Levy-Khinchin representation, we have (see e.g. Bertoin (1996); 
Sato (1999)) 

poo 

/ e-^'fni.)m^e-^''-^^\ (2.1) 
Jo 

where 

/•OO 

*o(s)= / {l-e-'-)7rn{du) (2.2) 
Jo 

is the Laplace symbol. First we provide an alternate proof of Theorem 3.1 of Meerschaert and 
Scheffler (2008). 

Theorem 2.1. For a strictly increasing subordinator Y{t) with density function p{x,t), the 
density function q{x, t) of the hitting time process W{t) — mi{x > : Y{x) > t} is given by 

q{x,t) = / TTvit - y, oo) fY{x){y)dy, t > 0, x > 0. (2.3) 

^0 

Proof. Since the sample paths of Y{t) are strictly increasing, we have P{W{t) < x) — 
P{Y{x) > t),x > Q,t > 0. Using the conditions q{x, 0) = = p(oo, x), we obtain 

rx poo 

/ q{u,t)du— / p{u,x)du, i > 0, x > 0. 
Jq Jt 
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This implies 



dxdt 



q{u, t)du 



dxdt 



p{u, x)du 



Taking LT on both sides and using Ctp{t,x) = e~^^'^^^\ we have 



Ct[q{x,t)] = -^y(s)e- 



d_ 

dx 



Ctpit, x) 



d_ 
dx 



(2.4) 



Using (12.2p and changing order of integration, we have 

e~*"7ry (u, oo)du = 



Tiyidy) 







TCyidy) )du 

y 



s 



(2.5) 



The inverse Laplace transform in f l2.4p can be written in convolution form of the inverse 
Laplace of \E'y(s)/s and e"^*^*^*-*. The result follows now by using the fact that e-^*^^') is the 
Laplace transform of the density function of the variable Y{x) and using (12. 5p . □ 

Let g{u,x) denote the density function of G{x) ~ IG(5x, 7) and let h{x,t) denote the 
density function of the hitting time process H{t). Here in the functions g and h the first 
and second variables are space and time variables respectively. Using (12. 4 p and the fact that 
Ctg(t,x) = e^''^(v^T^+2s-7)^ have the following result. 

Corollary 2.1. The LT of the density function of the hitting time process H{t) of the inverse 
Gaussian process G(t) is given by 

s 



Ct[h{x,t)] = ^(v/7T2;-7)e-'^^(v/^-7). 



(2.6) 



Remark 2.1. The Laplace- Laplace transform (LLT) of the function h defined by 



is given by 



h{u, s) 



h{u, s) 



e h{x,t)dt dx 



^(v/^^T27-7) 
sj M + (5(v/7^ + 2s-7)' 



(2.7) 



(2.8) 
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Theorem 2.2. The density function h{x,t) of H{t) can be put in the following integral form 
h{x,t) = -e^^''-'^ / — jsin{6x^/2^) + ^/2^cos{6x^/2^) dy. (2.9) 



Jo Z/ + 7V2 

Proof. The density function of H{t) is calculated by taking the Laplace inverse of £([^(^5^)] 5 
by using the Laplace inversion formula. Let Ct[h{x,t)] = F{s). Then 



1 



XQ+ioo 



h{x,t) = — e''F{s)ds. (2.10) 

■^'^^ Jxo-ioo 



Write 



F{s) = (5e^^" ^ v5^5e-^-y^ - l^-s^V^^j = Se^-"" {F^is) - -fF,{s)) (say). (2.11) 

For calculating integral in fl2.10p . we consider a closed key-hole contour ABCDEFA with a 
branch point at P= (—7^/2, 0). Here AB and EF are arcs of a circle of radius R with center 
at P, EC and DE are line segments parallel to a;-axis, CD is an arc 7,. of a circle of radius r 
with center at P and FA is the line segment from xq — iy to xq + iy with xq > 0. First, we 
calculate the above integral explicitly for 

FAs) = -e-^-V^\ (2.12) 
s 

By residue Theorem, 

/ e^^F,{s)ds + ^ [ e'^F,{s)ds + ^ [ e'^F,{s)ds + ^ [ e''F,{s)ds 
2111 Jab 27r? J 27rz Jcd 27r« J^e 

+ — : / e''Fi{s)ds + — : / e''Fi{s)ds = e-'^^^ (2.13) 



27ri Jef 27ri Jxo-iy 

which is the residue of Fi at s = 0. Since |-Fi(s)| < -, we have 

s 

lim / e'^Fi{s)ds= lim / e'^Fi{s)ds = 0, (2.14) 

(see Schiff (1988), Lemma 4.1, p. 154). Further, for the path CD, after putting s = —7^/2 + 
re*^, we have 



f e''F,{s)ds= [ —e-'^^^^'ds 
Jcd Jcd =^ 



^"""^ire'^de 

-72/2 + re*^ 

^0, (2.15) 
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as r — )■ 0, since the integrand is bounded. Along BC, put 7^ + 2s = 2ye™ so that + 2s 
i\f2/y and ds = —dy. We have 



._^2/2 St ^ — 

/ e^*Fi(s)c/s= / —e-^V^^+^^rfs 

r I e-'y-^'^^ydy (2.16) 

jR Ti^ + y 



' BC J-i?-7V2 

-472/2 / ^-ty-iSx 



Next along DE, take 7^ + 2s = 2ye this implies 1/7^ + 2s = —i^f2y and (is = — (iy. This 
implies 



/ e^*Fi(s)cis= / _e-^-V7^+2«(is 

= e-*^'/' ^77! e-*^+'^^^(i?/. (2.17) 



Using fl2l6l) and fl2Tfl) . we get 



/ e'^F^{s)ds + ^ [ e'^Fi{s)ds=''— f s%n{bx^y)dy. (2.18) 



Using fl27[3|) . (EH]), fl27[5|) and f l238|) with r ^ 0, i? ^ cx), we get 



— / e^*Fi(s)rfs = e-^^^ / -——sin{5x^)dy. (2.19) 

2vrzAo_ioo Jo 7/2 + 1/ 



Next consider 



F2(s) = -v/72 + 2se-'^"v/^. (2.20) 
s 

Using the same procedure as for -Fi(s), we get for F2(s) as 

The result follows now by using fl2A9D and f lCTl) with ( 12101) . □ 



. / e^*F2(s)rfs = 76-'^^ + / cos{5x^)dy. (2.21) 



Remark 2.2. (i) For the particular case 7 = and 6 = 1, 

V2 



\ 1 f"^ 

h(x,t) = — y-^/^e-'ycos{xJ2^)dy 
Jo 

/■oo 

/ e~*-^^ cos{V2xu)du 
Jo 



2V2 



'le-^y^\ t>0. 

TXt 



The last line follows by using the result /q°° e cos{au)du = \\f^(^ " ^"^^ (see e.g. Abramowitz 
and Stegun (1992)). 

(ii) It is obvious that B.{€) = |-B(t)|. However Hit) and |-B(t)| are not the same processes 
since the sample paths of Hii) are monotonically increasing however the sample paths of |-B(t)| 
are oscillatory. For a Levy process Z{t), E{Z{t)) = tE{Z{l)). Since tE{H{l)) = ^ 
yj^tji: = E{H{t)), we see that H{t) is not a Levy process. 

For g > 0, let Mg{t) = E{H{t)Y denote the g-th moment of the process H{t), which may be 
numerically evaluated for known t by using the density function h{x,t). However, an explicit 
expression can be obtained by using the LT of Mq{t). Let Mq(s) denotes the LT of Mq{t). As 
given in Veillette and Taqqu (2010a), we have 



where \I/g(s) = 5(a/7^ + 2s — 7) denotes the Laplace symbol of the process G{t). Now we have 
the following result for the mean and variance of the process H{t). Let Erf(;z) = ^ £ e"" du 
denotes the error function. 

Proposition 2.1. The mean function Mi[t) of H[t) is given by 

M.{t) = ^/^e--V2 + ^i?r/(7v^) + g (l + Erf[^^)) , 5, 7 > 0. (2.23) 
The second moment M2{t) of H{t) is given as follows; 

M,it) = ^ + ^ + e-^'*/2(7V^ + ^'V2t'/') + v^(2f t + ^H' - l)Erfi^^) 



52 (^2 p^2^ 

(2.24) 

For the particular case 7 = and 5 = 1, we have Mi(t) = 21/11 and M2{t) = 2t. 
Proof From dS^D, 

Mi(s) = , ^ . (2.25) 

5s(v/yT27-7) 

To invert this LT, we write 



Using the fact(see Roberts and Kaufman (1966), p. 210) 



^/s~+~a^ e 



+ yaErf(y^) (2.27) 



7 



and 



Ed{y)dy = zErf(z) + 



we have 



£-1 



1 



7 



dw + , 
'ixw a/2 Jo 



Erf(7v^/V2)dti; 



Also, £-i(l/s2) = t. Now the result follows for Mi(t) by using (|2:26|) with ([228]). 
Further, for M2{t) we have from fl2.22p 

1 



(2.28) 



Mo(s) 



525(^^2 + 25-7) 
1 

52 



Using similar steps as in (12.28^ . we obtain 



£-1 



+ 72/2 



^0 



7 

: dwdr H — ^ 

vrw ^2 



t rr 



JO 



Erf(7v^/V2)rfw;rfr 



273 ^/27r 



e-^ */2(7y2t + 73v^t3/2) ^ ^{2^H + - l)Erf(7yt72) 

(2.29) 



The result for M2{t) now follows using the fact that C^^{l/s^) = t^/2. 

For the particular case 7 = and 5 = 1, we have Mi{s) = s^^^^/v^ and JVhis) = 2/s^ which 
implies, Mi(t) = ^J^tpix and M^it) = 2t respectively. □ 

Next we have the following result regarding the asymptotic behavior of E{H(t)) and 
Var(if(t)). 



Corollary 2.2. As t ^ 00, we have 

E(H{t)) 

and 



^t, z/7>0 

1 /2t 
5' 



if 1 = 



1\2 



Further, as t ^ 0, 



Var{H{t)) ~ 



E{H{t)) ^\\ - and Var{H{t)) = o{t^'^) 

\ 71 



(2.30) 

(2.31) 
(2.32) 
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Proof. Using the fact that hm Erf(2) = 1, we have 



hm — — = hm 

t—>-oo t t—^oo 



1 e-'-i'/^ 



6V2^ Vi 257 1 



1 
6' 



Further, using the fact hmErff^;) = 0, we get 



hm — ^ = hm 



1 .e-*^^/2+ ^ ^ 



5V27r 



Erf (tv^) + (l + Erf (7 



+ ^ hm ^Erf ( -iJtTl 



1 /2 



1 1 

+ 



5V2^ 5V2^ 



The resuh for asymptotic behavior of variance of H{t) also foUows similarly. 



□ 



In the next result the LT of t) with respect to the space variable x is given in 
a closed form. However, this LT does not exist in a neighborhood of /i = and so it can't 
be used in finding the moments of H{t). Further, for 7 = 0, it can be used for finding the 
moments, since the LT exist for all /i > 0. 

Proposition 2.2. The LT (with respect to x variable) of the density function h{x,t) is 



C^[h{x,t)] 



-fix 



h{x, t)dx 



y2/i5e-^'/2 



TT 



-ty 



Proof. We have from f l2.9l) . 



Cx[h{x,t)] 



IT 



-ty 



2/ + 7V2L 



00 ^-ty 



Jo 

^b^f^Aj ^ yU — 57 



dy 



TT Jo y + 7V2 [ (/i - 57)2 + 252?/) 



+ ^2y 



(^-57)2 + 252^) 



TT 7o (l/ + 7V2)((/i-<57)' + 252y) 
The second last line follows by using the integrals 

h 



dy. 



e ""^ sinihx)dx = 
a2 + 6^ 



/■oo 

and / e~^''cos(bx)dx 
Jo 



a2 + b 



2,a>0,beM. (2.33) 



□ 
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Remark 2.3. For the case 6 = 1 and 7 = 0, 



n Jq /i2 + 2y 



''/'Erie 



(4) 



where Erfc{z) = 1 — Erf{z) 



Next, we look at the tail behavior of the density function h{x,t), which is useful for studying 
the infinite divisibility of the distribution function of H{t). 

Theorem 2.3. The tail probability P{H{t) > x) decay exponentially, as x —> oo. That is, 

P{H{t) >x) = 0(x-^e^^^-^'/^*), as X ^ oo. (2.34) 



Proof. We have from (12.91) . 



h{x,t) = -e^^^-^ 

71 



-ty 



s 

—e 

71 



1/ + 7V2 
27 



^sin{5x\/2y) + y'2ycos{6x\/2y) 



dy 



°° we-*" 



a;2 + 772 



sm 



{S\/2xuj)duj + 2V2 



00 ^2g-tc.2 



uj^ + lV2 



cos{S\/2xu)di 



(2.35) 



obtained by substituting u = ^/y. For a function q such that q{z) is real for real z and 
q{z) = Oiz-"^), r]>0, as Re{z) 00 (see e.g. 01ver(1974), p. 79) 



ze 



z^ + a- 



;q{z)dz = Trie °'^q{ia) + e(x), 



where 



ex 



ze 



z^ + 



q{z)dz 



(2.36) 



(2.37) 



with the line £: z = y + iw, w > a and —00 < y < 00. Taking q{z) = e and the 
corresponding ei(x), 



z"^ + 

Comparing the imaginary parts in both sides gives 



e-'^'dz = 7rze-""e*"' + ei(x). (2.38) 



^ ' e dz = ne + 
z^ + 

ZSin{xz) _j^2 71 _ax+ta^ 

^ , ' e dz = -e 



+ ei{x) 
+ eUx). 



(2.39) 
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Similarity, for q{z) = ze we get 



^ , / e *^d2; = --ae "V" +£2 0; . 2.40 

Also, ei(x) = e2{x) = 0{x^^e~^^^'^^), as x — )■ 00 (see Exercise 6.2 of Olver (1974)). Now using 
fl239|) and fICTD with (I235D . we have 

/i(x,t) = -(27 + 2y2)e^^"-^ei(x) 

= 0(x-^e^^"-"'/^*), asx->cx). (2.41) 

Now 



lim 



P(H(t)>x) , -h(x,t) 
lim 



< M, M > 0. 



[using L'Hopital's rule) 



Now the result follows. □ 

Remark 2.4. Since the tail probability P {H {t) > x) decay exponentially, moments of all order 
exist for h{x, t). 

A consequence of Theorem 2.3 is the following. 

Proposition 2.3. The distribution function of H{t) is not infinitely divisible. 
Proof. Since P{H{t) > x) < Mx~^e^'^^~^^^'^'^ for large x and some M > 0, we have 

2:^-00 y a; ma; / x^oo \ xmx I 

Hence, by Corollary 9.9 of Steutel and Van Harn (2004, p. 200), H{t) is Gaussian, which is a 
contradiction. □ 

Note that when 7 = 0, h{x, t) is a folded Gaussian density which is also not infinitely divisible, 
a known result (see e.g. p. 126 of Steutel and Van Harn (2004)). 

We next show that the hitting time of a stable process is also not infinitely divisible. 

Proposition 2.4. The tail probability of the hitting time process E{t) of the stable process 
D{t) has the following asymptotic behavior, namely, 

P(E(t)>a:) = 0(^^^^^i^x^-V^^^), as x^ 00. (2.42) 
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Proof. Let f{x, 1) denote the density function of a /3 (0 < /3 < 1) stable random variable 
D{1). Then 

m{x,t)^^x-^-^/^f{tx-^/^,l), x>0 (2.43) 

is the density function of the hitting time E{t) of the process D{t), defined by E{t) — mi{x > 
: D{x) > t} (see e.g. Meerschaert and Scheffler (2004)). Then, from (Uchaikin and Zolotarev 
(1999)) 

fix, 1) < K:r(i-/3/2)/(/3-i) exp(-(l - /3)(x/l3f/^^-^^), (2.44) 
for X sufficiently small. Hence, as x — > oo, 

mix, t) < K(te-i//5){i-/3/2)/{^-i) exp(-(l - p){tx-^/^/pf/^^-^^), (2.45) 
or in other words 

m{x, t)^0 (Mx(^-^/2)/(/5-i)-i-i/(/3) exp(-7Va;V(i-/3))) , (2.46) 
where M = (is:/^)t{i-^/2)/{^-i)+i and N ^ {I - p){t/ pf/^^'^l Thus, we obtain 



lim 

a;— >oo 



P{E{t) > x) 



lim 



(M(l - /3)/iV)a;(2-^/W-i)-V/3 exp(-iVa;V(i-/3)) 

-f{x,t) 



form 



(M(l - P)/N)£ (x(2-/3/2)/(/3-i)-i//3 exp(-ArxV(i-/5))) 

= 1. 

The result follows now. □ 
An application of the above result is the following. 

Proposition 2.5. The distribution function of E{t) is not infinitely divisible. 
Proof. Since 

P{E{t) >x)< ^(1^- /^) ^(2-/9/2)/(/9-i)-i//3 exp(-Ara;V(i-/^)) 
for large x, we have 

lnP{E{t) >x)< ln{M{l - j3)) - InN + ((2 - /3/2)/(/3 - 1) - l/j3)lnx - Nx^/^^'^K (2.47) 
Since, < ^ < 1, 

a:->oo \^ xlnx ) x-^oo y xlnx ) 

iVx'3/(i-/?) 

= lim ; = oo, 

x-^oo Inx 

and hence E{t) is not infinitely divisible. □ 
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The sample paths of IG process and its hitting time process H(t) are generated by using 
the algorithm given by Michael, Schucany and Hass (see Cont and Tankov (2004), p. 182 or 
Devroye (1986)). 
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Figure 1: Sample paths of IG process G{t) for 5 = 7 = 1 and the corresponding hitting time 
process H{t). 

In this section we obtain the partial differential equation satisfied by the density func- 
tion h{x, t). 

3 The PDE Connection 

Let for convenience we denote h{x, s) and h{u, t) as the LT and FT of h respectively. Now we 
have the following result. 

Theorem 3.1. The densities h{x,t) of the process H{t) solve 



(3,1) 



where 5o{t) is the Dirac's delta function. 



Proof. We have from ( 12. 6p 



h{x, s) = -(^7^ + 25 - ^)e-^"(V^-7). 



(3.2) 
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This implies 



—h{x, s) = -5(v/7' + 2s - ^)h{x, s), (3.3) 
ox 

which leads to 

—h{x, s) = 52(yyT27 - jYh{x, s). (3.4) 
Using dSJ]) and ((331) we get 

2^7— ) h{x, s) = 2S^{sh{x, s) - h{x, 0)) + 2S^h{x, 0). (3.5) 



^ dx^ dx ^ 

After inverting the above LT, we get 

d'^h ^ dh ^,dh 



a.^ 2.0^- = 2.^- + 2.^M..O)io«) (3,6) 

since £t((5o(t)) = 1. □ 

Remark 3.1. The density function g{x,t) of G(t) satisfies the following PDF (see Kumar et. 
al. (2011)) 

— - 26j- = 26 -. (3.7) 
Note that (13. ip has a closed connection with equation (13. 7p . 

The result in Theorem 3.1 can be generalized to the hitting times of tempered stable subordi- 
nators. It is well known that for a stable distribution with index < /3 < 1, the p-th moment 
is infinite for p > P, and its tail probability decays polynomially (see e.g. Samorodnitsky 
and Taqqu (2000)). To overcome this limitation, tempered stable (TS) distributions are in- 
troduced, which are infinitely divisible and hence correspond to a Levy process. Let f{x,t) 
denote the density function of a /3-stable process D(t) with index < /3 < 1 with LT 



oo 



f{x,t)dx = e-'' . (3.8) 

A tempered stable subordinator -D^(t) has a density f^{x,t) = e~'^^~^^'^^f{x,t), /i > 0. The 
Levy measure corresponding to a tempered stable process is given by (see e.g. Cont and 
Tankov (2004, p. 115)) 



-fix 



^D,Xx) = ^, OO, x>0. (3.9) 

POO 

Since / 7Cj;)^{x)dx = oo, Theorem 21.3 of Sato (1999)implies that the sample paths of D^{t) 
Jo ^ 

are strictly increasing, since jumping times are dense in (0, oo). Further 

/■oo 

/ e-'''f^{x,t)dx = e-'^^'+^^"'^"\ (3.10) 
Next we have the following result. 
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Proposition 3.1. The density function m^{x^t) of the hitting time process of D^{t) satisfies 
the following PDE. For (3 = \, n>2 

J2{-iy ('')f^^'-'^''^-^m^{x,t) = ^m,{x,t) + 6o{t)m,{x,0). (3.11) 
j=i V-^/ ^ 

Proof. We prove this result by induction. Using (12. 4p to get 

m^{x, s) = -{{s + fif - /)e-*«^+'^)'-'^'). (3.12) 



For n = 2, 



and 



dx'^ ' ^ 

Using (13131) and fl3AD . we get 

o,.i/2 d 



d 

— m^(x, s) = -{{s + nf/"^ - /i^/^)m^(x, s) (3.13) 
m,{x, s) = {{s + /x)i/2 - //^)H(^, (3-14) 



- 2/i^/^^ ) m^(x, s) = {sm^{x, s) - m^{x, 0)) + m^(x, 0). 



_ dx"^ dx 
Inverting the LT to get 



Similarly, for n = 3 



^ - 3//^^ + 3/i2/3- j m,(x, t) = (-1)^ (^"^/^(^' ^) + "^/^(^' 0) J ■ (3-16) 
The result also follows summarily for a general n. □ 

Further we have following results on fractional partial differential equation. Let /dt^ denote 
the Caputo fractional derivative defined for < /3 < 1 by 

1 /■* f'(r) 

w^m - I (3.17) 

where /' denotes the derivative of /. 

Proposition 3.2. The density function h{x,t) of H{t) satisfy the following pseudo-fractional 
partial differential equation 

d^h{x, t) + 5 (^^72 + 2dt - 7) h{x, t) = 0. (3.18) 

For 5 = 1 and 7 = 0, the density function h{x,t) satisfies 

^Hx,t) + V2^h{x,t) + \j^h{x,^) = 0. (3.19) 
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Proof. Using (jM]), the LIT of h{x,t) satisfies 

Lh{u, s) + (5(^7^ + 2s - 'j)h{u, s) = s~'^S{^/'j'^ + 2s - 7). 



MA 



Since h{0, s) = s ^5(a/7^ + 2s — 7), we have 

uh{u, s) - h{0, s) + (5(^7^ + 2s - j)h{u, s) = 0. 
Invert the LT with respect to variable u to get 

dMx, s) + 5(^7^ + 2s - -f)h{x, s) = 0. 
Now invert the LT with respect to variable s to obtain 

dM^, t) + S [Vl^ + 25* - 7) h{x, t) = 0. 
Further, for 5 = 1 and 7 = 0, we have from (I3.20p 

uh{u, s) + \/2sh{u, s) 

which implies 

uh{u, s) — h{0, s) + \/2sh{u, s) = 0. 
Invert the above LT to obtain 

d r 

Thus 



— hiu, s) + \/2sh(u, s) = 0. 
ox 



d ~ ( ~ 

— s) + ^2 f y/sh{u, s) 



Inverting the above LT and using (see e.g. Podlubny (1999), p. 106) 



0. 



(3.20) 



(3.21) 



s^Ctifix, t)] - s^~'f{x, 0), < /3 < 1 

the result follows. 

Remark 3.2. The density function of G{t) has FT (see e.g. Applebaum (2009), p. 54) 
Then 



d_ 

di' 



(u,t) = -(5(a/72 + 2m - -f)g{u,t), 



after inverting the FT 



dtg{x, t) = -6{^/^^ + 2^^ - -f)g{x, t), 



(3.22) 

(3.23) 
(3.24) 

(3.25) 
(3.26) 

(3.27) 
□ 

(3.28) 

(3.29) 
(3.30) 
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and which has a close connection with f l3.22p . In particular for 6=1 and 'y = 0, we have 
from fl3.30l) the LT with respect to the variable x is 

gM = e-'^'. (3.31) 
Using similar lines as in Proposition 3.1 and the fact that g{0,t) = 0, we have 

d r- d^/^ 

Q^9{x,t) + V2^g{x,t) = 0, (3.32) 
and which is again closely related with (13.321) . 

Subordinated processes have interesting connections to PDE (see e.g. Allouba (2002); Allouba 
and Zheng (2001); Baeumer et. al. (2009)). We consider here the subordinated process 
X{t) = B{H{t)), where H{t) is the hitting time of the inverse Gaussian process and B{t) is 
the standard Brownian motion. Let u{x,t) denotes the density function of the process X{t). 
Let hx denotes the partial derivative of h with respect to x. 

Proposition 3.3. The density function h{x,t) satisfy the following properties. 

lim h{x, t) = h{0, t) = Se-^'/^ - je^''^^Erfc (j^^ j • (3-33) 

Further, 

lim hJx, t) = hJO, t) = 26-fh(0, t) (3.34) 
and, as X ^ oo, hn,{x,t) = 0{Nxe^^''~'-''^ ^'^^) , for some N > 0. 

Proof Note h{x,t) = _e^T^-TV2 / I(t;x,y)dy, where I{t;x,y) follows from Then 

Jo 

\Iit'.x,y)\ = —^—^['ysin{6x^/2y) + >/2ycos{6x^/2y)] 



< 



I/ + 7V2 



^Jly^i^sinibx^lA) + B) , (where sin{e) = (3.35) 



which is integrable. By dominated convergence theorem, we have 
lim h{x, t) = -e-^'/2 / 0^ y^dy 
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Note that I(t;x,y) and 
d 



-ty 



I{t;x,y) = - 
ox y + Y/^ 



6'^ \/ 2ycos{6x y 2y) — 26'^sin{6x\/2y) 



(3.36) 



are continuous in x and y. following the steps for (13.351) . 



Q^nt;x,y) 



2^bJye-^y , 
< — - — - — - — (independent of x) 

2l/ + 7 



(3.37) 



and / — — 7: — dy < 00. Hence, / —I(t;x,y)dy is uniformly convergent and so we 

Jo "^y + T Jo ox 



can differentiate under the integral sign. Thus, 



5 



and 



TT 



57 / I{t;x,y)dy + 
2?/ + 7 



d 



I{t;x,y)dy 



dx 
^dy = 257/1(0, t). 



(3.38) 



Using similar argument as in Theorem 2.3, we can show that, as a; — )■ 00, 

/ —I(t; X, y)dy = 0(Cxe''^^-^'/^*), for some C > 0. 
Jo 

Thus, h^{x,t) = 0(A^xe''^'^-'^'/^*), using fICTD and flCTjl along with I KM . 
Theorem 3.2. The densities u{x,t) of the process X{t) satisfies the following PDE 

2<52-m(x, t) = + 5-f—y{x, t) + ft{x, 0)/i(0, t) - 25\(x, Q)5,{t). (3.40) 



(3.39) 
□ 



Proof. Since the inner and outer process are independent, the density function of X{t) is 

u{x,t) = I w{x,r)h{r,t) dr, (3-41) 
Jo 

where w{x,t) is the density of the Brownian motion B{t), and h{r,t) is the density of the 
process H(t). Since the FT w{u,t) = exp(— ^tw^), we have 

dw 
'dt 

and hence 

dw 
'dt 



(ivS^w 



1 d^w 

2'dx^ 



(3.42) 



(3.43) 
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is the governing equation for B{t), the usual diffusion equation. Now write 

d f°° d 

26 —u{x, t)= w{x,r)2S^—h{r,t)dr 



(9^ d \ 

w{x,r)(—-— — 26'j—)h{r,t)dr — 26'^6o{t) I w{x,r)h{r,0)dr 







d 



d 



dr 



w{x,r)—h{r^t) ~ —w{x,r)h{r,t) — 26'yw{x,r)h{r,t) 



=0 dr 



r=0 



r=0 



roo q2 roo q 

+ / -TT-^w{x.ir)hir,t)dr + 25'^ / —w{x,r)hir,t)dr — 25'^u{x,{i)5Q{t) 
Jo or Jo or 

-w{x, 0)/i^(0, t) + wt{x, 0)/i(0, t) + 257w(x, 0)/i(0, t) 

-j^ i-co qA poo q2 

H — / -——w{x,r)hir,t)dr + 5'^ I —^w(x,r)h(r,t)dr — 26^u(x,0)6o(t) 
4 Jq dx'^ Jq dx^ 

-w{x, 0)h^{0, t) + wt{x, 0)h{0, t) + 26-fw{x, 0)h{0, t) 

+ + b^^^u{x,t) - 252w(x,0)(5o(t) (using ^M) 

and hence the resuh follows. □ 
Remark 3.3. For the 'particular case 5=1 and 7 = 0, we have the following fractional PDF. 

Ql/2 roo Ql/2 



+ (57— t) + 0)/i(0, t) + 2(5^m(x, 0)5o(t) (using (fSJl) 



^1/2 ^00 ^1/2 

V 2^^m(x, ^) = y ^(2^' '")'^2^^/i(r, t)rfr 



00 



9 / 2 

w(x,r)—h(r,t)dr — \ — / w(x,r)h(r,0)dr 
•9'" V TTt Jo 

= w{x,0)h{0,t) + - / —rw{x,r)h{r,t)dr-\ —u{x,0) 
2 Jq dx^ V vrt 

1 92 AY 
= w{x, 0)/i(0, t) + ^Q^uix, t) - y — m(x, 0). 

For detailed and interesting discussion on Brownian time Brownian motion B{\B{t)\), see 
Alloba and Zheng (2001) and Beaumer et al. (2009). 
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